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Abstract 

We evaluate the pion mass in finite volume to two loops within Chiral Perturbation Theory. 
\ The results are compared with a recently proposed extension of the asymptotic formula of 

Liischer. We find that contributions, which were neglected in the latter, are numerically 
very small at the two-loop level and conclude that for M-j^L > 2, L > 2fm the finite volume 
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o 

effects in the meson sector are analytically well under control 

o 
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1 Introduction 

Numerical simulations in lattice QCD are bound to rather small volumes. When 
one determines the hadron spectrum and other low energy observables one has to 
understand and properly account for the volume dependence in order to correctly 
interpret the numerical data. Analytical methods which allow one to predict the 
size of the finite volume effects are particularly useful in this respect. In the case 
of hadron masses, there are two different methods to analytically evaluate the 
finite volume effects: the asymptotic formula derived by Liischer and chiral 
perturbation theory (ChPT) in finite volume jSEllll- Liischer's formula relates 
the volume dependence of the mass of a hadron to an integral over the vr-hadron 
forward scattering amplitude. Knowledge of the latter scattering amplitude im- 
mediately translates into an estimate of finite volume effects for the hadron mass. 
Since the integral is dominated by the low-energy region, one can rely on the chiral 
representation of the scattering amplitude in order to numerically evaluate the in- 
tegral. In this manner one makes use of ChPT only in infinite volume and obtains 
an estimate only of the leading exponential term in the finite volume dependence, 
the term of the order exp(— MttL), where L is the box size. Alternatively one can 
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perform the calculation of the hadron mass in ChPT in finite volume and obtains 
an estimate also of the terms which are exponentially subleading. 

A detailed study and comparison of the two approaches has been performed 
in 1^1 E] for the case of the pion mass. This case is particularly interesting both 
because it is a simple hadron to be studied on the lattice and because the vrvr 
scattering amplitude, which is needed in the Liischer formula, is now known to 
next-to-next-to- leading order [71 |H1 13 El in the chiral expansion. It was therefore 
possible to study numerically how well the leading exponential term dominates 
the series, and how fast the chiral expansion (of the leading exponential term) 
converges. The somewhat surprising result was that the leading term in both series 
receives large corrections both from subleading exponentials as well as from the 
next-to-leading order in the chiral expansion. The next-to-next-to-leading order 
chiral correction, on the other hand, was found to be rather small, which indicates 
that - at this level ~ the chiral series has started converging well. 

The fact that the Liischer formula appeared to give the numerically dominating 
term only for volumes so large that the finite volume correction itself has become 
negligible appeared to be more worrisome. At first sight this might have lead one 
to conclude that in cases of practical interest one could not rely on this extremely 
convenient formula in order to evaluate reliably the finite volume effects. We 
believe, however, that such a negative conclusion is unjustified and have proposed a 
resummation of the Liischer formula |lH ll2j which retains all the convenience of the 
original one but does not suffer from the same large corrections. This resummation 
can be understood in simple terms: Liischer has shown that the leading exponential 
correction comes from a radiative correction in which the emitted virtual light 
particle (the pion, in QCD) goes around the world once (thanks to the periodic 
boundary conditions) before being reabsorbed. The resummed formula takes into 
account all other possible ways which the pion has to go around the world (go 
along the diagonal, or go around more than once, etc.). 

As was shown in ^21 this resummed formula exactly reproduces the one-loop 
ChPT calculation of the pion mass in finite volume if one inserts in the integral the 
leading chiral representation of the vrvr scattering amplitude. Inserting the next- 
to-leading chiral representation of the scattering amplitude one reproduces the full 
two-loop calculation of the pion mass in ChPT in finite volume up to corrections 
of order exp(— ML) with M > {y/S + 1)/\/2Mt,-. Despite this improvement of the 
algebraic accuracy of the resummed formula with respect to the Ltischer's one, 
it is essential to show that the resummation is numerically effective and that the 
corrections to it are small. In order to make this check we have now performed 
a full two loop calculation of the pion mass in finite volume in ChPT. As we 
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will demonstrate in this article, the corrections which are not captured by the 
resununed asymptotic formula are negligibly small for M^L > 2. The results of 
the present two-loop calculation were anticipated in jl3lll4| . Here, we give further 
details about the calculation and the results. 

One of the reasons for performing this further investigation of finite volume 
effects for hadron masses is that one may view Liischer's formula (or its resummed 
version) as a way to determine on the lattice a scattering amplitude, cf. ^j. While 
it is true that the scattering amplitude is seen here in an exponentially suppressed 
effect, it is also true that the direct calculation of a scattering amplitude on the 
lattice is much more difficult than that of a mass. This indirect method to extract a 
scattering amplitude may in some cases turn out to be more practical, and we find 
it worthwhile to discuss its theoretical feasibility. Indeed if the part of the finite 
volume corrections which is related to the scattering amplitude is not strongly 
dominating with respect to the rest, this is not a viable method to determine the 
scattering amplitude. The conclusion reached in this paper is therefore encouraging 
in this respect and can be used to estimate in which region of the {Mt^,L) plane 
the finite volume correction to a hadron mass is given to a good approximation by 
the integral containing the scattering amplitude. 

The two-loop calculation is interesting in its own right, also from a technical 
viewpoint. To date, a number of finite volume calculations have been performed 
at one-loop order ^1 El E[ El El 1201 , but as far as we know the only finite 
volume two-loop calculation which has been performed until now is for the quark 
condensate jUj, which is a much simpler calculation. In the related context of 
finite temperature field theory there are a few examples of calculations beyond 
one loop HSIHH. 

We wish to briefiy mention related work. Most notably, the asymptotic formula 
may also be applied to the nucleon mass [15] , see Koma and Koma (2^1 as well as 
jl4j for recent work. Braun, Pirner and Klein have evaluated the volume depen- 
dence of the pion mass based on a quark-meson model |27j . In the framework of 
a lattice regularised ChPT finite volume effects have been addressed by Borasoy 
et al. in Ref. |2H]. 

The outline of the article is as follows. In sect.[21we set the notation and remind 
of the basic assumptions for an application of ChPT in finite volume. Sect.|31is 
devoted to outline the calculation and state the main results. In sect. I^we show 
the explicit expressions which have been used for the numerical analysis in sect. [SI 
We conclude with a summary. 
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2 Preliminaries 



In this section we shall set the notation and the basic definitions for the two-loop 
calculation. 

2.1 ChPT in finite and in infinite volume 

Chiral Perturbation Theory (ChPT) is the effective theory for QCD at low energies. 
If we first restrict ourselves to the infinite volume case, the effective Lagrangian of 
QCD for two light fiavours at low energies consists of an infinite number of terms 

El, 

C^s = C2 + U + Cq + ... . (1) 

As we wish to calculate the pion mass, an on-shell quantity, external fields can be 
dropped in >Ceff- We work in the isospin symmetry limit rriu = rrid in Euclidean 
space-time, and for the choice of the pion fields we use the non-linear sigma model 
parameterization. We have 

4 



^2 = — {u^u^ - x+) , (2) 



with 







X = 25ml , m=^{mu + md) , (3) 

with = U. The symbol (A) denotes the trace of the two-by-two matrix A. The 
Lagrangian £4 can be written as [7j, 

4 

1=1 

where 

1 1 1 i 

Pi = -^{ut^Uf^f , P2 = --{Uf,u^f , P3 = - J^ix+f , Pi = -^{u^X-ti) , (5) 

with 

= u^d^xu^ - ud^X^u . (6) 
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The ellipsis in eq. (0]) denotes terms that do not contribute to the pion mass. 
The low energy constants ii are divergent and remove the ultraviolet divergences 
generated by one-loop graphs from £2- 

The complete effective Lagrangian Cq with its divergence structure at d = 4 
has been constructed in [^H EHI- As will be discussed in sect. |31 terms from Cq 
merely renormalize the pion mass in infinite volume and do not contribute to finite 
volume corrections which we are interested in. Thus, we refrain from showing it 
here. Given the effective Lagrangian and the parameterization for the pion fields, it 
is straightforward to calculate the pion mass to two-loops. We refer to [S], where 
one also finds a detailed discussion of the renormalization procedure [two-loop 
diagrams in ChPT are discussed in jHlj]. 

The effective framework is still appropriate, when the system is enclosed by a 
large box of size V = L^. We refer to the literature for the foundations and 
a recent review jllj. Here, we only recall a few fundamental results which guided 
the present calculation: the volume has to be large enough, such that ChPT can 
give reliable results, 2Ft^L 2> 1. The value of the parameter M^L determines the 
power counting for the perturbative calculation: if M-,^L ^ 1 one is in the "p- 
regime" in which 1/L counts as a small quantity of order M^^. If Mj^L < 1 one is 
in the "e-regime" and is a quantity of order M,r- In both cases the effective 
Lagrangian is the same as in the infinite volume. In this article we only consider 
the "p-regime" , where the system is distorted mildly and the only change brought 
about by the finite volume is a modification of the pion propagator due to the 
periodic boundary conditions of the pion fields^ 



with Gq{x) the propagator in infinite volume. 
2.2 Basic definitions 

In Euclidean space-time the propagator is defined through the connected correla- 
tion function 



where a,h stand for isospin indices of the pion fields and the subscript L in the 
correlation function denotes that it is evaluated in finite volume. 

^Throughout we denote by the volume the three-dimensional volume V — Li^ , whereas the 
time direction is not compactified. 




(7) 




(8) 
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We also have 



where the momenta p can only have discrete values 

27]- 

p=— n, nGZ^, (10) 

and = IBih is the tree-level pion mass in infinite volume. The pion mass in 
finite volume M-^l is now defined by the pole equation 

G(pl)-i=0, for pL = {iM^L,0) . (11) 



3 Outline of calculation and statement of results 

For a large volume, the finite size effects are expected to be small, such that the 
pole equation can be solved perturbatively. We outline the calculation and proceed 
with the main results. More details about the calculation are relegated to later 
sections or the appendix. 

3.1 One— loop result 

Since the effective Lagrangian remains unchanged when going to the finite volume, 
we can immediately write down the Feynman diagrams which contribute to the 
self-energy at two-loops, see fig. ^ the only difference with respect to an infinite 
volume calculation is that the propagators need to be periodified, cf. eq.Q- At 
leading order, the graphs 1(a) and 1(e) need to be evaluated and the self-energy 
admits the form 

Sl(/,p) = ^G(0)(-^M,^-P^)-2^3^ + o(^) , (12) 

with G(0) the value of the finite volume propagator at the origin. It contains a 
logarithmic divergence due to the contribution from the term n = in eq.((7|. In 
dimensional regularization, 

Go(0)^pip/<iV^^^^M-. (13) 
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(a) (b) (c) (d) 



(e) (f) (g) (h) 



Figure 1: Self-energy graphs to two-loops in ChPT. A spline corresponds to a 
periodified propagator, whereas those without correspond to an infinite volume 
propagator, cf. eq.Q. Normal vertices come from L21 squared vertices from 
and the circle-crossed from L^. 

The remaining terms with n 7^ are finite and may be expressed in terms of a 
kinematical function^ (7i(M^,0,L), 

G(0) = Go(0) + 5i(m2,0,L), 

For a derivation of ea. H14() . we refer to [32] • In app.^we provide a different deriva- 
tion, based on a contour integration analysis. The pole of the Gamma function 
in ea.H13|) in four dimensions is absorbed in a renormalization of the low-energy 
constant £3. One readily verifies that inserting ea.lfTHl into ea.lfT^ yields for the 
leading finite volume shift [2j 



l + ^5i(M^,0,L) + O(-l 



(15) 



The separation of the cut-off and the volume-dependence is as expected: finite 
volume corrections do not generate new uv-divergences. At leading order the finite 
volume corrections could be isolated immediately. This will not be the case at the 
two-loop level. The graph (d) in fig. ^ does not factorize in pure one-loop integrals 
and further steps need to be performed. We refer to app. ^ for further details. 



^The second argument of g\ is the temperature, which we keep zero. Notation as in Ref.[5| 
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3.2 Minimal set of periodified propagators 

In the evaluation of the Feynman diagrams, it is mandatory to make use of the 
fohowing result. In order to evaluate a graph consisting of L loops, one needs to 
consider only a certain set of L propagators as the periodified, finite- volume ones — 
the others can be taken as infinite-volume propagators. In the following we briefly 
discuss this statement (see also p.lSff in Ref.P). Consider an arbitrary self-energy 
graph with L loops, / internal lines ^ and V vertices. Since the number of loops 
is the number of independent integrations over momenta, we have 

L = I-V + l. (16) 

For every line ^ the propagator is an infinite sum of terms characterized by an 
integer vector n(£). The loop graph itself is an infinite sum of terms which can 
be identified by a set of integer vectors {n(£)} assigned to all the internal lines. 
As remarked by Liischer this structure can be viewed as a gauge field on the 
self-energy graph. Gauge transformations are defined as'^ 

n'(^) = n(£)+A(/(^))-A(i(£)), (17) 

where A(t;) is some field of integer vectors and i{t) {f{i)) is the initial (final) 
vertex of the line i. One verifies immediately that two contributions of a Feynman 
graph which differ only by a gauge transformation yield the same mathematical 
expression. In the calculation of the loop graph what matters is the sum over 
representatives of gauge equivalence classes. A convenient representative can be 
found, e.g., by adjusting A{i{£)) and A{f{i)) iteratively, such that n{i) = for as 
many internal lines i as possible. This can be achieved for V — 1 lines. Therefore, 
there remain I — V + 1 internal lines where the periodified propagator has to be 
inserted - which coincides with the number of loops of the graph. This shall be 
our minimal set of periodified propagators. In fig. ^ we have attached a spline to 
a periodified propagator, whereas lines without a spline correspond to an infinite 
volume propagator. 

3.3 Two— loop result 

It is convenient to split the sum over the equivalence classes into three parts'^: 



S(2) 



n(^) = V ^ (pure gauge), 

n(^) = V ^ except for one line £ (simple gauge), (18) 
n(£) = y i except for two lines £1,^2- 



^We borrow the notation from Ref.0. 

^Note the sUght difference in the definition of simple fields with respect to T*. We do not 



require \n(£)\ 
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The pion mass in finite volume to two-loops then admits the form 
= Ml- S(^) - , 
Ml = m2-S(°) , (19) 
where, using = M^L, we get 

E« = Ip + I, + 0{e) , (20) 

oo 



n=l 



-oo 



oo oo , — ; TT- 

^0 = -3^1:^ ^S^^^^d.c[^.(M + i.)l . (22) 



with m{n) = number of integer vectors z with = n. Ip denotes the pole and 
Ic the cut contribution, whose meaning and precise definitions will be explained 
below. The expression for is more cumbersome: 



= mI(, 



^9i{Kf - ^A^5i(A^)^5i(A^) + A 



+ Oie) , (23) 



and we merely note that it can be split into products of one-loop contributions 
(terms with and a pure two-loop part, denoted by A. We have introduced the 
abbreviation 

A detailed derivation of these results will be given in the subsequent sections and 
we confine ourselves to a few comments at this stage: In Ip one recovers the 
asymptotic formula of Liischer, if one restricts the sum to the first addendum. Its 
extension to the present form has already been suggested in and was applied 
in ^21- The function .7v(z^) denotes the isospin zero tttt scattering amplitude in the 
forward kinematics. It contains cuts due to the two-pion intermediate state. In the 
derivation of the asymptotic formula, Liischer consistently dropped contributions 
arising from the cuts, since in a large volume expansion they are beyond the order of 
accuracy he aimed at. Here, we take them into account up to two-loops, relying on 
the chiral representation of vrvr scattering amplitude. The outcome of the analysis 
is summarized in Ic in which the same symbol this time with two arguments 
appears: with the latter we mean the scattering amplitude of two on-shell pions 
into two off-shell pions in the forward kinematics configuration. The contributions 



9 



from the cuts can still be written as integrals over the vrvr scattering amplitude with 
an exponential weight, as in the asymptotic formula. Notice however, that this 
formula explicitly relies on the chiral representation of the scattering amplitude 
and is only valid up to two-loops, whereas the Liischer formula holds to all orders. 

Contributions from two pion propagators in finite volume are ultimately cap- 
tured in E(2), being expressed in terms of a dimensionless function giiX-K) and a 
numerically small correction A arising from graphed). Both are explicitly given 
in aPD. \K\ and eq. respectively. 

The pure gauge contributions are not volume dependent and merely renormal- 
ize the pion mass, cf. eq. H19() . A detailed discussion of this calculation can be 
found in IH] . with which we agree - this was a useful, nontrivial check for our 
calculation. 

3.4 Self-energy to first order: E^^) 

The simple fields can be summed up in closed form and may be represented by a 
skeleton diagram, see fig. |2h) 

S(^) = 2/ (^E"^(")«'''^'^^o('z')r..(p,g,-p,-g) , (25) 

with p = {iM-jr, 0) and r7r^(p, q, —p, —q) the 4-point function of vrvr scattering in 
the forward scattering kinematics. Note that the 4-point function in eq. H25() is an 
off-shell amplitude and that we evaluate it at p and not at pL- The difference is 
taken into account in T,^'^\ The representation of S*^^) in eq. (|25jl has already been 
used by Liischer and eventually led him to the asymptotic formula He then 
discussed the contribution of the pole of the propagator Go{q'^) that one meets at 

qi = i^M^ + ql+qi, q± = {q2,q3). (26) 

Above this pole the singularities come from the cuts of the propagator and the 
4-point function T.,rnip, q, —p, —q)- These start from 

s = -{p + qf >iM^ , u= -ip - qf > 4M^ , -q^ > 9M^ . (27) 

Liischer showed that the pole contribution is dominating with respect to those 
coming from the cuts and neglected the latter. In fact, his discussion involved 
no further assumptions about the 4-point function and remains true at every 
order of the perturbative expansion. Since our goal is to test the asymptotic 
formula beyond the leading exponentials, we wish to work out the impact of the 
contributions that were dropped by Liischer. Doing this at the two-loop level is 
rather straightforward, as we will now show. 
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a) 



b) 



Figure 2: Skeleton diagrams representing a) eq. H25|) and b) eq. 1)64(1 . The blob in 
o) stands for the 4~point function of vrvr scattering in infinite volume and in b) for a 
subtracted 6-point function of Trvrvr scattering in infinite volume. The double-line 
with the spline is a finite volume propagator with the physical pion mass M^. 



Up to O(p^), the four-point function can be decomposed into a combination 
of functions which have either a singularity in s or in u. Since S*-^^ is symmetric 
in s and u, we may write r,r7r(p, q, —p, —q) in terms of a function r^^(s,p(7), which 
depends only on the variables s and pq, and which has cuts at s > 4M^, 

r7r7r(P,g,-P, -g) = r^7r(s,M) + ^-^^ • (28) 

We write qi = x + iy and for the domain s G M, s > 4M^ we have 

-x^ -ql-ql> 2,Ml , xy = -qoM^ . (29) 

In particular we observe that in the complex qi plane the cut starts above the 
pole of eq. (|26() . as illustrated in fig. |31 Performing now the contour integration 
in the upper half plane we obtain two terms: the first one, Ip, comes from the 
pole and is simply its residuum. The other comes from the integral along the new 
integration path. The latter contribution vanishes as we push the integration lines 
to infinity, except for the path which goes around the cut, to be denoted by Ic in 
the following, 

=Ip + Ic. 

As was shown by Liischer the former can be simplified considerably and can 
be brought in the usual form of the asymptotic formula 



n=l 



— oo 



11 





Imqj 




c i 








Reqj 



Figure 3: Integration contour in the complex qi plane with the pole from the pion 
propagator and the branch cut from the vrvr scattering amplitude. 

with T-,^{iy) the tttt forward scattering amphtude. Restricting the sum to the first 
term we recover Liischer's formula pOj. For Ic we find 

I^ = _Y,m{n) j^^^ J dqr^^^disc[fUs,pq)]+0{F-^), (30) 

[s>4A/2] 

where disc[r7r7r(s,M)] denotes the discontinuity of F^^ along the cut. It is now 
convenient to shift the integration path in qq from Im(qo) = to Im(go) = — iM^r. 
Along this path we have 

90 = 90 - iM^ , s = -Qq - qf - qj_, qo (31) 

and the integration over qi = x + iy falls onto the imaginary axis, 



X = 0, y>yQ = ^AM^ + ql + ql , 

oo oo _ 



yo 



n=l 

Next, we change the integration variable from qi to s 



4M2 



n=l 



and make use of 

(27r)2 2(/i2 + g2)i/2^ 47r^L ' 

to carry through the integration over q± and to end up with eq. ((2H). 
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3.5 Self-energy to second order: S*^^) 

The term S^^^ is more complicated to manipulate, although it is suggestive to 
think that it may be related to an integral over the Svr Sir amplitude. A similar 
question has been asked in the framework of finite temperature QCD 23 . We 
discuss the situation in finite volume in appendix ^ Since these considerations 
have not lead us to a nice and compact representation for E^^^ we have written it 
as in H23() where we only have split the factorizable two-loop contributions from 
the rest. The latter, which we have denoted by A we have evaluated numerically. 

We conclude this section with a remark concerning the large L behaviour of 
the two-loop results. In the large volume limit, the contributions from eq. H21jl 
behave according to 

lim /p ~ ^ e~^^ . (35) 

Similarly, we may evaluate the large volume behaviour of the terms occurring in 
eqs. (|23l22|) . It turns out that these, besides being exponentially suppressed and 
behaving at least as exp(— aAvr), with a = (\/3 + l)/v^, are also suppressed by 
a power of Ajr- Unfortunately this suppression is not particularly strong, being 
of I/i/A^ with respect to eq. (|35|) . While we can conclude that the resummed 
Liischer formula is dominating in comparison to other two-loop diagrams in a 
1/V%r expansion, i.e. 



M^L = - , (36) 

we are convinced that the most important test of the resummed formula is the 
numerical one, which we will discuss in the following. 



4 Summary of analytical results 



In this section we shall give the analytical results in explicit form, i.e. insert the 
chiral representation of the amplitudes which appear in the formulae given in the 
previous section and express our results in terms of a few basic integrals. In ea. 1)191 - 
we have split the finite size effects of the pion mass into Ti^^^ = Ip + Ic and 
T,^'^\ For the former two we find 



n=l 



m{n) 1 
n Ayr 



r(2) 
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Ir 



n=l 



n Ayr 



(37) 



where ^ is the chiral expansion parameter defined in ea. (|24l) and the expressions 
-^mI' -^mI have already been given in Ref. [01^1, and we reproduce them here for 
completeness: 

r(2) _ dO 



^(4) 



-B 



(38) 



55 ,7j 8^ 5^ 

h 4^1 + -£2 £3 

18 3 2 



2£a 



112 8^ 
~9~ ~ 3 



32 



, 13 16 1 40 2 



16 

y 



Re 
Im 



B^{x) = 2K2{x)/x 



k even 
k odd 



where 

and the integrals i?o ^^'^ defined as 

1 

dy e" V "(^+^') 5(2 + 2iy) for 

where g{x) is related to the standard one-loop function J{xM^) through 

g{x) = 16TT^JixM^) , 
and J{xM^) given in eq. (|56j) . For instance for x < 0, we have 

with 



g{x) = crlog^— ^ + 2 



CT+1 



a = x/l - 4/x , 



(39) 



(40) 



(41) 



(42) 



and elsewhere defined through analytic continuation. In Ref. also the next-to- 
next-to-leading order term in Ip, namely ^'^I^^ has been given and we will use it 
in our numerical analysis. The coefficient Ic can be expressed as a combination of 
basic integrals 

ic = - (112C0'2 + 37Ci'0 - 40Ci'2 - 4C2'0) , 
3 



dy I d~s g2 + 4^2 f 1 



1/2 



(43) 



Notice that the latter expression is obtained from the chiral representation of the 
off-shell vrvr scattering amplitude for forward kinematics. This is parameterization 
dependent, and the result given here is obtained for the parameterization discussed 
in sect. 12. ll The dependence on the parameterization must cancel in the full result. 
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The self-energy to second order has already been introduced in eq.(|231) and repro- 
duced here for convenience 



with 



A 



fievr 



2^2 



(44) 



where H,H^ and H^j^i, are related to the sunset-type integrals of fig. E^). We 
have not been able to find a compact representation for these integrals and only 
elaborate on their numerical analysis in app.^in some detail. 



5 Numerics 
5.1 Setup 

The numerical analysis is performed in line with the setup of Ref. \12\ . The quantity 
of interest is 

Rm. = , (45) 

whose quark mass dependence shall be evaluated numerically for different sizes L. 
The parameters of Rm^, are (see eqns.(^H21l and eqns.(j3Z[S])) the pion mass 
and the pion decay constant F.^ in infinite volume as well as (implicitly in /|^^^^) 
the SU(2) low energy constants. The quark mass dependence of the pion decay 
constant may be taken into account by expressing as a function of the pion mass 
(see, e.g. [H]). Regarding the low energy constants, we use the ones determined 
in [^1 ^] which are the same as in our previous finite size studies jSl El ■ 



5.2 Results 

We plot our results for -Rm^ in figs.EJ and El both for L = 2,3,4fm as a function 
of and for = 100, 300, 500MeV as a function of L. We show the one-loop 
result (LO) as well as the two-loop result (NLO). These shall be compared with the 
resummed asymptotic formula with LO/NLO/NNLO input for the vrvr scattering 
amplitude. Note that the one-loop result and the resummed asymptotic formula 
to LO coincide. The best estimate for Rm^, is finally obtained by adding to the 
asymptotic pure three-loop contribution the two-loop result (NNLO asympt. + 
full NLO). At NLO, the finite size effects contain low energy constants, see eg. 
diagram /) and g) in fig. ^ leading to a non-negligible error band. 
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M^(GeV) 

Figure 4: Rm^ = M^l/M^ - 1 vs. for L = 2,3,4fm. The result of the 

resummed asymptotic Liischer formula ()21|) with LO/NLO/NNLO chiral input 
(with attribute "asympt." for NLO and NNLO in legenda) is compared to the 
one-loop (LO) and two-loop (NLO) result. The best estimate for Rm^ is obtained 
by adding the pure three loop contribution from the asymptotic formula to the 
two-loop result (NNLO asympt. + full NLO). The error band comes from the 
uncertainties in the low energy constants and is only shown for the best estimate. 
In the region above the Mt^L = 2 line, one is not safely in the p-regime and our 
results should not be trusted. 




L(fm) 



Figure 5: Rm^ = M^l/M^ - 1 vs. L for = 100, 300, SOOMeV. The rest of the 
legenda is as in fig. ^ 



L = 2 fm 


= 140 MeV 
M^L ~ 1.4 


= 200 MeV 
M^L ~ 2.0 


= 250 MeV 
M^L ~ 2.5 






LO 


0.0463 


0.0199 


0.0105 




Ip 


ANLO 


0.0291 


0.0149 


0.0088 




Ic 


ANLO 


0.0005 


0.0001 


0.0000 




S(2) 


ANLO 


0.0206 


0.0038 


0.0011 




Ip 


ANNLO 


0.0076 


0.0053 


0.0035 


Total 


0.1041(55) 


0.0440(50) 


0.0239(42) 



Table 1: -Rm^ = ^I-kl/J^-k — 1 for a 2 fm volume and pion masses M,r = 140MeV, 
200MeV, 250MeV. We show the relative numerical impact of the LO, the pure 
NLO and the pure NNLO contribution for Rm^- In the second column we give the 
source of the effect. The fifth line, e.g. contains the contribution to Rmt^ of the 
dispersive terms Ic alone, once kinematical prefactors are added accordingly. The 
numerical results for the Ip contributions, to LO, NLO and NNLO have already 
been given in [HI Il2] . 

We take up a point which was already observed in [3, namely the large con- 
tributions when going from LO to NLO in the asymptotic formula (dotted to 
thin~dash~dotted) . Compared to this gap, the additional contributions from the 
full two-loop result (thick-dash-dotted) are very small. The two-loop and the 
NLO result from the asymptotic formula only drift away, when we go beyond the 
region where the p-regime can be safely applied. In tab. 1 we wish to underline 
these statements with a numerical example: we show the relative numerical im- 
pact of the LO, the pure NLO and the pure NNLO contribution for RM-n- Iii ths 
second column we give the source of the effect. The fifth line, e.g. contains only 
the contribution of the dispersive terms Ic to -Ra/^i once the kinematical prefactors 
are properly accounted for. We observe that Ic is strongly suppressed, irrespective 
of the value M^^L. Consider the column with Mt^L ~ 1.4. Although the bulk of the 
subleading effects is still due to the asymptotic contributions, the terms of S^^^ 
play a significant role and can not be neglected. This behaviour was expected, 
since with M^-L ~ 1.4 we might have already crossed the border of the p-regime. 
As we increase this parameter to M-^-L ~ 2, the asymptotic regime begins to set 
in. The contributions from the resummed asymptotic formula are now dominating 
with respect to those from S^^^. The numerical results for M-,^L ~ 2.5 confirm 
this trend. The fact that even the additional NNLO asymptotic terms (a par- 
tial three-loop result) are larger than the NLO non-asymptotic contributions is 
in nice agreement with the analytical expectation found in eq. (|H6j) . However, we 
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find that the suppression seen in the numbers is stronger than what one could have 
expected on the basis of the latter analytical argument. FinaUy, the discussion of 
the subleading effects allows us to give a reliable estimate for the lower bound of 
M-j^L for the p-regime, 

Mtt-L > 2 : lower bound for p-regime . (46) 

6 Summary 

i) We have evaluated the finite volume corrections for the pion mass to two- 
loops within the framework of chiral perturbation theory (ChPT) in the 
p-regime (M^L > 1,L > 2fm, < SOOMeV). 

a) We have compared the two-loop result with the resummed version of the 
asymptotic formula of Liischer. We have found that whenever the effects 
are calculated for masses and volumes such that M^^L ^ 1, such that one 
is safely within the p-regime of ChPT, the contributions which are not in- 
cluded in the resummed asymptotic formula are very small. The result gives 
us confidence in the claim that the resummed asymptotic formula is a conve- 
nient and efficient way to reliably calculate finite volume effects for hadronic 
masses jT2]. 

in) The derivation of the asymptotic formula for decay constants follows 
closely the original one for the masses - the only new, subtle point, 
concerns the amplitude which appears in the integrand of the asymptotic 
formula, which has a pole in the integration region which first needs to be 
subtracted (see for details). This is, however, a technical point - the 
physics of the finite volume corrections for masses and decay constants ap- 
pears to be rather similar. In view of this we believe that the present results 
speak also in favour of the resummed asymptotic formula for decay constants 
|12j . A check of this claim could be obtained by evaluating the pion decay 
constant to two loops. 

iv) The two-loop calculation performed here allows us to better estimate the 
region of validity of ChPT in the p-regime. Our conclusion is that in the 
case of the pion mass it is necessary to have M-,^L > 2. Again, this serves as 
a guideline also for decay constants and masses of other hadrons. 
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A Finite volume integrals 

In this appendix we give further details on the finite volume integrals which oc- 
curred in the two-loop calculation. Throughout we applied dimensional regular- 
ization, as is common in ChPT. Since the finite box breaks Lorentz invariance, 
tensor simplifications have to be performed with care. Consider e.g. the integral 

WV = y(^L Mf^2(P<l) ' (47) 

which arises when evaluating diagram /) in fig. ^ Remember p = {iM-,r,0). The 
prime in the sum denotes that the term with n = is excluded. The ansatz 

A^u = 5^uA (48) 

with A a scalar integral which is valid in infinite volume, leads to an incorrect result 
in finite volume. A direct calculation of the integral with the help of eq. H61() yields 

^^.^^ V = -^^S^ E ^^./2(A.y^) . (49) 

In four dimensions the result agrees with the derivation outlined in the next section, 
where one proceeds along the same lines as in sect. \\\A\ 



A.l Tadpole 

In fig. n the finite volume corrections of the diagrams (a)-(c),(/) and {g) factorize 
into one-loop integrals which are of the generic form 

with k an integer positive number and V{pq,q'^) a polynomial of pq and g^. It 
suffices to discuss the case for /c = 1, since k = 2,3,... are obtained through 
appropriate derivatives with respect to M^. The only singularity of the integrand 



19 



is the pole of the propagator. Therefore, the contour integration analysis applied 
in sect. 13.41 yields the result 



- / dy e- V"(i+s/')^- ViiM^y, -Ml) 

n J 



^ 87r2A. 

n=l * ^ 

In particular, in the text we have used the dimensionless function 51 (At^), 
167r2 



(51) 



gkiMl,0,L) 



Ml 

E' 



■gi{Ml,0,L) , 



d'^q 



^iqnL 



{2TT)d {Ml + (?2)fc 



(52) 



evaluated at d = 4, and where gi{Ml,0, L) was introduced a long time ago by 
Gasser and Leutwyler |2l Oj , see also eq. (fH]) . 

A. 2 Sunset 

The only real two-loop diagram is the sunset graph fig. ^) , and we will comment 
on it in some detail. It is convenient to split the finite volume integrals as in 
ea.(|18|) only after the tensor simplifications. As alluded in the beginning of the 
appendix, tensor simplifications in finite volume have to be performed with care. 
Even though one can not rely on Lorentz invariance, the sunset tensor integrals 
may still be reduced to the structures 



{n,n^,n^,}= / d'^xe^^^Gixf {i,ia^,-a^ajG(x) 



(53) 



with G{x) from eq. ((TJ. A rather direct way to perform these steps is to work in 
coordinate space and to make use of partial integrations, i.e. 



d'^xe'P''d^,G{x)d^G{x)Gia 



I d'^xe'P'^Gix'^ 



Pf,iduG{x) + d^duG{x) 



■ (54) 



Notice that the same identities in momentum space are derived with the help of 
translational invariance which is still respected in finite volume due to the periodic 
boundary conditions. In the following, we will only elaborate on the scalar integral 
Ti. in more detail. We expand the integral in terms of number of finite volume 
propagators as motivated in sect. 13. 3[ 

n = + 3W(^) + 7^(2) , (55) 
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where the first (second) addend corresponds to the pure (simple) gauge fields 
contribution. For TC^'^^ we refer to [S^. Further, 

W«=^m(n) J j^,jjr^,J[{p-qf]+9i{MlO,L)m , (56) 



n=l 

with 

„2\ 



(27r)'^ [M2 + £2] [M2 + (A; - £)2] 

The first term of eq. (|56j) is finite, the second carries an uv-divergence which is 
absorbed by a counterterm. This shows that although finite volume effects do 
not generate new uv-divergences, they still appear at intermediate steps of the 
calculation. It is a thorough check on our calculation that these non-analytic 
divergences cancel. Finally, 

^(2)^ V f "'^"^ "^'^'^ 1 ^.o^ 

^ / (27r)M2vr)MM2+g2][M2 + A2][M2 + (^_g_fc)2] " ^ ) 

In the text, we used its dimensionless version 

^ ' " " J (2-)^ (2.)^ [1 + [1 + ^2] [1 + (p - ,~ - ~ky] ' ^''^ 

which are uv-finite and do not need to be renormalized. We only have to find a 
convenient representation for the numerical analysis. We shall restrict ourselves 
to the scalar integral in the following. Introducing a Feynman parameter by com- 
bining the second and third denominator, we find 



H 



We use Schwinger's trick for both remaining denominators 

oo 

1 



dae-"(i+^ ) . (61) 



The integrals over k and q are then of the Gaussian type and can be performed 
analytically. We are then left with three integrations over a rather lengthy ex- 
pression which shall not be written down here. Despite their unhandy form, the 
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integrations may still safely be performed numerically. The accuracy of the deter- 
mination of eq. H6U|) is not limited by the integration routine, but by the rather 
slow convergence of the sum in n and r for moderate X-j^. Consequently, the eval- 
uation of the sunset integrals going into A is restricted to three significant digits. 
(The last digit given for in tab. 1 is not significant.) Note that the uncertainty 
of the H-type integrals is not a serious matter. Firstly, it could be lowered by 
brute force and secondly it only plays a (minor) role, in case when the p-regime 
can not be safely applied anymore. 

B Self— energy to second order: E^^^ 

In this section we ask ourselves whether the self-energy to second order can be 
represented in a similar compact form as in the case of the self-energy to first 
order in eq. As will be discussed, it is indeed possible to relate the self- 

energy to second order to a Svr scattering amplitude in the forward scattering 
kinematics, as illustrated in fig. |2t). This has already been observed by Schenk in 
a related context ||22j. He investigated the dynamics of pions in a cold heat bath 
of temperature T and examined the effective mass of the pion M^(r) within the 
framework of ChPT. The close relation between finite temperature field theory and 
finite volume effects becomes apparent in the imaginary time formalism, where one 
treats the inverse temperature as a finite extension in the imaginary time direction. 
The expansion in terms of the number of finite volume propagators is then in 
one-to-one correspondence with the expansion in terms of the number of finite 
temperature propagators. In the following, Schenk's approach shall be adapted to 
the finite volume scenario. We first establish the relation between S^^^ and the 
three-to-three particle scattering amplitude and proceed with various remarks. 

Consider the 6-point function in d dimensions in the non-linear sigma model 
parameterization eq. © in the forward kinematics 

3 

^ y dxi . . . dx5e-i^'(---^)-i'=(-^--=)-'''(-3-x6)(o|T^i^^-^(^b^V,i^^^^(^^jo^ 

a,b=l 

" (M2 + p2)2(^2 + P)2(Af2 + ^2)2 ^--(^'' ^' ' (^2) 

with Z the wave function renormalization constant and ipx = fix). The amplitude 
T-n-TTirip, k, q) contains a pole at p'^ = —M^ which needs to be subtracted 

TnnniP, k, q) = f^j,n{p, k, q) + ^^f ' ^' '^^ . (63) 
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bbbbb b a b 

a) b) c) d) e) 



Figure 6: Decomposition of the Stt-Stt amplitude in terms of Iparticle irreducible 
parts. The characters 1, a and b on the external legs denote isospin indices. 

The self-energy to second order can then be written in terms of the subtracted 
amplitude TTTTnrip, k, q) 

8 J (2^)4 (2^)4 (M2 + + A;2) ^ 

where the prime restricts the sum to integer vectors n and r obeying n 7^ 7^ r, 
and for diagram c) in fig. El in addition n 7^ r. The latter restriction avoids double 
counting as this term is already accounted for in a simple gauge field of T,^^\ 

i) The reason for the subtraction is easily accounted for. In fig.|Hlwe decompose 
the Svr-Svr scattering amplitude in terms of 1-particle irreducible parts. The 
subtraction removes the contribution from the diagram in fig. ^) which 
does not correspond to a 1-particle irreducible self-energy diagram, once 
the external legs are appropriately closed. 

A thorough discussion on the physical interpretation of the pole term can be 
found in ^5] . 

ii) Notice that eq. defines a subtracted off-shell amplitude which depends 
on the regularization scheme as well as on the parameterization of the pion 
fields. At the order we are working, the regularization dependence is not 
an issue, since the subtracted amplitude T.,^Trn{p,k,q) is only needed at tree 
level. However, the dependence on the parameterization of the pion fields is 
of concern. While the momentum integrations in eq. (|64j) for the diagrams 
fig.int) andinii) put the momenta k and q on-shell and the ambiguity due to 
the parameterization therefore drops out, this does not happen in the case 
of diagram Et) . 

Note that the same parameterization ambiguity already occurred in Ic in 
the dispersive analysis of S^^^ (cf. eq. pT|) where 7^ —M^). In order to 
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understand the close relation between these two terms, we first note that only 
diagram fig.^d) contributed to Ic- Further, the simple gauge field of fig-^^d) 
with n 7^ for one propagator can immediately be written as a contribution 
with two finite volume propagators: one periodifies a second propagator, 
however with the same n 7^ as already for the first one. Since Mj^l does 
not depend on the parameterization of the pion fields, the ambiguities of 
the two terms have to cancel each other. To explicitly show this is however 
nontrivial. 

In summary, while our representations for Ic in eq. (|22|) . resp. eq. (|37|) and 
for S(2) in eq. (jHU, resp. eq. do depend on the off-shell dependence of 
the scattering amplitudes, the sum in eq. (fT^ does not. 

Hi) Even though the self-energy to second order can be expressed in a compact 
form, further simplifications (similar to those performed in sect. 13.41 seem 
not to be straightforward. Instead, for the (numerical) integrations we had 
to discuss the various terms contributing to S^^^ one by one. A general 
discussion of finite volume integrals has been provided in app. After 
treating the finite volume integrals, we end up with the basic functions given 
in eq. (|^. 
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